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X ■ Abstract 

We compute the next-to-next-to-leading order (NNLO) contributions to the three splitting functions 
governing the evolution of unpolarized non-singlet combinations of quark densities in perturbative 
QCD. Our results agree with all partial results available in the literature. We find that the correct 
leading logarithmic (LL) predictions for small momentum fractions x do not provide a good esti- 
mate of the respective complete results. A new, unpredicted LL contribution is found for the colour 
factor d"^'^dabc entering at three loops for the first time. We investigate the size of the corrections 
and the stability of the NNLO evolution under variation of the renormalization scale. Except for 
very small x the corrections are found to be rather small even for large values of the strong coupling 
constant, in principle facilitating a perturbative evolution into the sub-GeV regime. 



1 Introduction 



Parton distributions form indispensable ingredients for the analysis of all hard-scattering processes 
involving initial- state hadrons. The dependence of these quantities on the fraction x of the hadron 
momentum carried by the quark or gluon cannot be calculated in perturbation theory. However, 
the scale-dependence (evolution) of the parton distributions can be derived from first principles 
in terms of an expansion in powers of the strong coupling constant ttg. The corresponding nth- 
order coefficients governing the evolution are referred to as the n-loop anomalous dimensions 
or splitting functions. Parton densities evolved by including the terms up to order a"+^ in this 
expansion constitute, together with the corresponding results for the partonic cross sections for the 
observable under consideration, the N"LO (leading-order, next-to-leading-order, next-to-next-to- 
leading-order, etc.) approximation of perturbative QCD. 

Presently the next-to-leading-order is the standard approximation for most important processes. 
The corresponding one- and two-loop splitting functions have been known for a long time HI |2l 
|3l|4l|5l|6l|71|8l|9l[l0l[TTl. The NNLO corrections need to be included, however, in order to arrive 
at quantitatively reliable predictions for hard processes at present and future high-energy collid- 
ers. These corrections are so far known only for structure functions in deep-inelastic scattering 
lfT2lfT3lfT4lfT5l and for Drell-Yan lepton-pair and gauge-boson production in proton-(anti-)proton 
collisions lfT6l [TTl [TSl fT9ll and the related cross sections for Higgs production in the heavy-top- 
quark approximation ifTTlllOl 121112211 . Work on NNLO cross sections for jet production is under 
way and expected to yield results in the near future, see Ref. ll23l and references therein. For 
the corresponding three-loop splitting functions, on the other hand, only partial results have been 
obtained up to now, most notably the lowest six/seven (even or odd) integer-A'' Mellin moments 

These Mellin moments already provide a rather accurate description of the splitting functions 
at large momentum fractions x Il25ll27ll28ll29ll . Their much-debated behaviour at small values of 
X, on the other hand, can only be determined by a full calculation. As we will demonstrate below 
for the non-singlet cases, this statement holds despite the existence of resummation predictions 
for the leading small-jc logarithms ll30l l3T1l . since -a- the correctly predicted logarithms do not 
dominate the three-loop splitting functions at any practically relevant value of x and -b- a term of 
the same size occurs with a new colour factor at third order which could not have been predicted 
from lower-order results, analogous to the situation for the four- loop (3-function of QCD ll32l . 

In this article we present the (unpolarized) flavour non-singlet (ns) splitting functions at the 
third order in perturbative QCD. The corresponding flavour singlet results will appear in a forth- 
coming publication O^ . The present article is organized as follows: In section 2 we set up our 
notations for the three independent third-order splitting functions and briefly discuss the method 
of our calculation. The Mellin-A'' space results are written down in section 3 together with their 
explicit large-A'^ limit which is relevant for the soft-gluon threshold resummation ll34l l35l l36ll at 
next-to-next-to leading logarithmic accuracy ll37ll . A surprising relation is found between the lead- 
ing large-A/^ term at two loops and the subleading {\nN)/N contribution at third order. In section 
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4 we present the exact resuhs as well as compact parametrizations for the x- space splitting func- 
tions and study their behaviour at small x. The numerical implications of these results for the 
scale dependence of the non-singlet quark distributions are illustrated in section 5. Except for very 
small values of x, the perturbation series appears to be well-behaved even down to sub-GeV scales 
where the initial distributions have been studied using non-perturbative methods for example in 
Refs. Il38ll39ll40ll4ni42ll43l . Finally we briefly summarize our findings in section 6. 



2 Notations and method 

We start by setting up our notations for the non-singlet combinations of parton distributions and the 
splitting functions governing their evolution. The number distributions of quarks and antiquarks in 
a hadron are denoted by and respectively, where x represents the fraction of the 

hadron momentum carried by the parton and stand for the factorization scale. There is no need 
to introduce a renormalization scale Hr different from fJ-f dX this point. The subscript i indicates the 
flavour of the (anti-)quark, with / = l,...,n^ for flavours of light quarks. 

The general structure of the (anti-)quark (anti-)quark splitting functions, constrained by charge 
conjugation invariance and flavour symmetry, is given by 

P — P - — 8 ; + P^ 

■^q/qt = -^q/q*: = ^ikPqq+Pqq- (2.1) 

In the expansion in powers of tts the flavour-diagonal ('valence') quantity starts at first order, 
while /'qq and the flavour- independent ('sea') contributions P^^ and P^^ are of order a^. A non- 
vanishing difference P^^ — P^^ occurs for the first time at the third order. 

This general structure leads to three independently evolving types of non-singlet distributions: 
The evolution of the flavour asymmetries 

= ± - i^k ± ^k) (2-2) 
and of linear combinations thereof, hereafter generically denoted by q^^, is governed by 

^ ns ^ qq -^^ qq • 

The sum of the valence distributions of all flavours, 

"/ 

^ns= K^.-^.), (2.4) 

evolves with 

pv _ pV _ pV I /pS _ ps \ _ p I ps (1 
ns ^ qq qq ^ f\ qq qq^ ~ -* ns ^ -* ns • ) 

The first moments of P^^ and P^^ vanish, since the first moments of the distributions q^ and 
reflect conserved additive quantum numbers. 
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We expand the splitting functions in powers of tZs = OCg/ (4;r), i.e. the evolution equations for 
q^^{x,ijj), z = ±, V, are written as 

^ = j P W'(^) ® ^^^(^,^2) (2.6) 

where ® represents the standard Mellin convolution. 

Our calculation is preformed in Mellin-A^ space, i.e., we compute the non-singlet anomalous 
dimensions Jm'iN) which are related by the Mellin transformation 

Yif(A^) = -/'rfxx^-ipif(;c) (2.7) 
Jo 

to the splitting functions discussed above. The relative sign is the standard convention. Note that 
in the older literature an additional factor of two is often included in Eq. (12.71 ). 

The calculation follows the methods of Refs. EH EH SH. We employ the optical 
theorem and the operator product expansion to calculate Mellin moments of the deep-inelastic 
structure functions. Since we treat the Mellin moment A^^ as an analytical parameter, we cannot 
apply the techniques of Refs. lEHEHEill' where the Mincer program ||35l|371 was used as the 
tool to solve the integrals. Instead, the introduction of new techniques was necessary, and various 
aspects of those have already been discussed in Refs. P31H51B91l5m . Here we briefly summarize 
our approach, focussing on some parts which have not been presented yet. It should be emphasized 
that we have at our disposal a very powerful check on all our derivations and calculations by 
letting, at any point, N be some positive integer value. Then we can resort to the approach of 
Refs. lEHEIlEil and, with the help of the Mincer program, the checking of all programs greatly 
simplifies. 

We start by constructing the diagrams for the forward Compton reactions 

quark (P) + vector (2) — > quark (P) + vector (Q) , (2.8) 

which contribute to the non-singlet structure functions Fj, and F3 of deep-inelastic scattering. 
The A^-th Mellin moment is given by the A'^-th derivative with respect to the quark momentum P 
at P = 0. The diagrams are generated automatically with the diagram generator Qgraf fSTl and 
for all symbolic manipulations we use the latest version of FORM ll5^ 15^ . The calculation is 
performed in dimensional regularization ll5?l 1551 15^ 15711 with D = 4 — 2e. The unrenormalized 
results in Mellin space are formulae in terms of the invariants determined by the colour group ll5Sll . 
harmonic sums E [71 |52l 1^ EUl and the values ^3, ^4, ^5 of the Riemann ^-function. In physics 
results the terms with ^4 cancel in A'^-space. With the help of an inverse Mellin transformation 
the results can be transformed to harmonic poly logarithms Il62ll63ll64l in Bjorken-;c space. Details 
have been discussed in Refs. P31l651l . The renormalization is carried out in the MS-scheme ll66ll67l 
as described in Ref. Il24ll25ll26t . 
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The complete non- singlet contributions to the structure functions can be obtained from three 
Lorentz projections of the amplitude for the process (I2.8I) . that is with g^'^, P^'P^ and with £^2/'v ^ 
ga(3/jvp^gp For the projection with g-"^ and P^'P^ one has two vector- like couplings, whereas 
for the projection with £^2a'v qj^^ j^^g jj^g product of a vector and an axial-vector coupling. The 
axial nature leads to the need for additional renormalizations with Z^, the axial renormalization, 
and with Z5, the finite renormalization due to the treatment of the 75. This is all described in the 
literature f6^. For the anomalous dimensions we need only the divergent parts of the g^^ and 
£^2a'v projections, but just as for the fixed moments we can also obtain the finite pieces which lead 
to the coefficient functions in N^LO. The determination of the latter for F2 and requires also 
the computation of the pt^P^ projection which is still in progress. The results for the three-loop 
coefficient functions will thus be presented in a future publication ll69ll . 

To solve the integrals we apply the following strategy lH51l^ . We set up a hierarchy of classes 
among all diagrams depending on the topology, for instance ladder, Benz or non-planar. Within a 
certain topology, we define a sub-hierarchy depending on the number of P-dependent propagators. 
We define basic building blocks (BBB's) as diagrams of a given topology in which the quark 
momentum P flows only through a single line in the diagram, while composite building blocks 
(CBB's) denote all diagrams with more than one P-dependent propagator. We determine reduction 
schemes that map the CBB's of a given topology class to the BBB's of the same topology class or 
to simpler CBB topologies. Subsequently, we use reduction identities that express the BBB's of a 
given topology class in terms of BBB's of simpler topologies. 

This procedure has been discussed to some extent in Refs. lH51 H^ . It exploits various cat- 
egories of relations between the integrals which can be derived as follows. For a generic loop 
integral depending on external momenta P and Q, the first category are integration-by-parts iden- 
tities lElEDl, 

/n^''^"a^^i'x(--) = 0- (2.9) 

These give a number of nontrivial relations by making various choices for the pi and pj from the 
loop momenta. Additionally pj can be equal to P or Q. The second category is based on scaling 
arguments P3I in Mellin space. They involve applying one of the operators 

,, d ,, d ,, d 

Qf"- — , P^^ — , P^"^ — (2.10) 

both inside the integral and to the integrated result. The scaling in Mellin space tells us the effect 
of these operators on the integrated result, while inside the integral we just work out the derivative. 
These relations naturally involve polynomials linear in A'^. The fourth operator of this kind, 

2^3^, (2.11) 

cannot be used naively in this context, because it does not commute with the limit P P 0. More 
care is needed in this case and we will come back to this shortly. 
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A third category of relations is obtained along the lines of the Passarino-Veltman decomposi- 
tion into form factors [71 1. In Mellin space we write 



/ ]\d''pnP'^ X (...) = Q^Q + PHp , (2.12) 

where Iq and Ip are the two form factors. By contracting Eq. (12.121) either with or P^, the Iq 
and Ip are determined in terms of a number of integrals. Next, by taking the derivative with respect 
to Q^, the relevant identities can be obtained. Because the momentum p, can be any of the loop 
momenta, Eq. (12.121) gives us as many relations as there are loops. Again, in Mellin space, these 
relations contain polynomials linear in A^. 

The fourth and the fifth category of relations are new. Together with the form factor relations 
from Eq. (12.121) they were crucial in setting up the reduction scheme for the three-loop topologies. 
They are based on operators that do not commute with the limit P P ^ 0. In the fourth category, 
one considers the dimensionless operators 

PQ „ 9 

01 = — -Q^^, (2.13) 

QQ dP^' 

d a 

02 = PQ^^, (2.14) 

03 = - — (2.15) 

Q Q dP^' dP^' 

Each individual operator Oi does not commute with the limit P P ^0, but certain linear combina- 
tions of the Oj do. However, one has to extend the ansatz based on scaling arguments in A''-space. 
Specifically, one has for the A'^-th moment of an integral I{N) 

w ^ (||)"(e.arcr.(f|)""^(a.ercf..... (2.16) 

where the C^^ and are dimensionless functions of A^, and a adjusts the mass dimensions. 
The novel feature is here the term ' proportional to P ■ P, which one may call higher twist. In 
contrast, for the relations based on Eq. (12.101) it was sufficient to restrict the ansatz to 

Applying the differential operators Oj in Eqs. (12.131 ) - (12.151 ) to the ansatz (12.161 ). one finds that 
the combinations 



2(a+l-A^)Oi-02 , (2A^-4 + D)C>i-03 (2.17) 

do commute with the limit P P ^ 0. That is to say, any dependence on the higher twist term 
Cf^^ vanishes in this limit and one is left with only contributions from C^^ . Eq. (12.171 ) adds two 
more relations, which in Mellin space contain quadratic polynomials in A'^ due to the differential 
operators of second order. We have checked that differential operators of yet a higher order in P 
and Q do not add any new information. 
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Finally, the fifth category of relations again uses the form factor approach of Eq. (I2.12I) . How- 
ever, now we do not take the derivative with respect to but with respect to P^. Some extra 
book-keeping is needed here, since one has to take along terms proportional to P P. Let us write 
Eq. (I2l2b as 

pfl = Q^lQ+P^'Ip. (2.18) 
Taking the derivative of Eq. (12.181) with respect to P^ in //-space one finds 

^^f^ = Q'^lQ + iD + N-l)Ip. (2.19) 



Solving Eq. ( 12.181) for Iq and Ip as above, however keeping all terms P P, substituting into 
Eq. (12.191 ) and finally taking the limit P P 0, we find 

5 u P-Pi „ d , .Q PiP Q-Q QP Pi 

.pi^I = — ^Q^'^/-K(D-fA^-2) ^ ^' ^ — —I. (2.20) 



Again, as the momentum pi can be any of the loop momenta, Eq. (12.201) gives us as many relations 
with polynomials linear in N as there are loops. 

Taken together, the reductions of category one to five suffice to obtain a complete reduction 
scheme. In particular, the reduction equations of category two to five involve explicitly the param- 
eter A'' of the Mellin moment. They give rise to difference equations in A'' for an integral I{N), 

ao{N)I{N)+ai{N)I{N-l) + ...+a,n{N)I{N-m) = G{N) , (2.21) 

in which the function G refers to a combination of integrals of simpler topologies. Zeroth order 
equations are of course trivial, although sometimes the function G can contain thousands of terms. 
First order difference equations can be solved analytically in a closed form, introducing one sum. 
Higher order difference equations on the other hand can be solved constructively, sometimes with 
considerable effort, by making an ansatz for the solution in terms of harmonic sums. For the 
present calculation we had to go up to fourth order for certain types of integrals. 

Due to the difference equations, which have to be solved in a successive way, a strict hierarchy 
for topology classes is introduced in the reduction scheme. For a given integral /, a difference 
equation as in Eq. (12.211 ). with some (often lenghty) function G expressed in terms of harmonic 
sums, can be solved in terms of harmonic sums again. Subsequently, the result for / can be part 
of the inhomogenous term in a difference equation for another, more complicated integral. This 
requires the tabulation of a large number CBB and BBB integrals, because each integral is typically 
used many times, thus it saves computer time and disk space. Only this tabulation, which required 
the addition of features to FORM ll5^ . renders the calculation feasible with current computing 
resources. For the complete project, including Refs. 0^1^ . we have collected tablebases with 
more than 100.000 integrals and a total size of tables of more than 3 GByte. 
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3 Results in Mellin space 



Here we present the anomalous dimensions yts^{N) in the MS-scheme up to the third order in the 
running coupling constant ttg, expanded in powers of as/ (47r). These quantities can be expressed 
in terms of harmonic sums |l6l Ul EH |60||. Following the notation of Ref. ll59ll . these sums are 
recursively defined by 

S±m{M) = £ ^ (3.1) 



i=l 



and 



S±mi,m2,...,mi,{M) = ^ — Sm2,...,mi,{l } 



(3.2) 



i=l 



The sum of the absolute values of the indices m^t defines the weight of the harmonic sum. In the 
n-loop anomalous dimensions written down below one encounters sums up to weight 2n — l. 

In order to arrive at a reasonably compact representation of our results, we employ the abbre- 
viation Sffi = S0i{N) in what follows, together with the notation 



N±5r„ = S,^{N± 1) , N±i5,^ = Sa{N±i) 



(3.3) 



for arguments shifted by ±1 or a larger integer i. In this notation the well-known one-loop (LO) 
anomalous dimension EJEIl reads 



y'SiN) = Q(2(N-+N+)5i-3) , 
and the corresponding two second-order (NLO) non- singlet quantities ||4l|6| are given by 

fi 



(3.4) 



Y^;^"(A^) = 4QC^ (2N+53 - ^ - 25_3 - y^i + (N_ + N+) + 2Si,_2 - 



+ 4C^n^(j^ + ^5i-(N-+N- 



+ 4C/,2(^45_3+2Si+2S2 



S2 + 253 



(N_+N- 



5i+45i,-2 + 251,2 + 252,1 +53 
Si- S3 



8 



2(N_+N+-2)5i 



(3.5) 



(3.6) 



The three-loop (NNLO, N^LO) contribution to the anomalous dimension 7+ (A^^) corresponding 
to the upper sign in Eq. (12.31) reads 



Yi?+(A^) = 16QC^n/(|C3-^ + y5_3-y53 + |5i,_2-|5_4 + 25i,i-^52 



10 



25 



257 ^ 2 ^ 

1237 11 
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52 1 53 1 54 

. ' 3 ' 3 



(Nh 
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23 
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(N_+N- 



^1,1 



16 2 1111 
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1 



1 



1883 



11 3 

+ ^"^-3,1 + 2*^-2 - 65-2^3 - 25-2-3 + 35-2,-2 - 45-2-2,1 + 85-2,1,-2 

13 



3 ' 2 
10Si,-3 3 

+ (N- +N+ - 2) [351C3 + 1 15i,i -45i,i,-2 
91 



■5i 



z 

16 5 1 176 13 

- ^^1,-2 + 125l-2,1 + 45l3 - 452,-2 - -^S^ + -^S^ + —52 + y 53 
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' ' " - +(N-+N - 



9131 19 
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- 27^1 + 9^2 + (N- + N+) [-51 - -52 + -53J ) + 16C/Q (- 
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+ (N- + N+-2) [65iC3-y5i+355ia-125i,i,_2 + 5i,2 + 1052,-2 + 52,1 +252,2 -253,1 



355 



+ (N_+N+) 
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(3.7) 



The third-order result for the anomalous dimension Jj^^iN) corresponding to the lower sign in 
Eq. (12.31 ) is given by 



yI? (N) = Yi?+(A^) + 16QQ fQ-^)((N_+N+-2) i:^5i + 1251^3 +25i, -2 
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1255 
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853 -2 + 2654 



(3.8) 



Finally the quantity y^f^iN) corresponding to the last term in Eq. (12.51 ) starts at three loops with 
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+ -(N+-N+2) 



54+52,-2—53,1 



■(N-2+N+2) 



5l,-3 — 5i,-2,l — 5i,i,-2 



+ (N--1) 



1 



454 + -55 



+ (N--N- 



1 1 

x52 —3 + -zSl —2 ■ 

2 ' 2 ' 



109 41 67 

^52 --52,1 + -53 



1 113 

^53,1 -52,1,-2 + ^52,3 + 2*^3,-2 - ^54,i 



y5i-^5i,-3+25i,-2,i-5i,i,-2) . (3.9) 



Eqs. (13.71) - (13.91) represent new results of this article, with the exception of the (identical) 



parts of Eqs. (13.71) and (13.81) which have been obtained by Gracey in Ref. I.72J and of the 

9 



contribution linear in «^ in Eq. (13.71) which we have published before All our results agree 
with the fixed moments determined before using the Mincer program ll46l l47ll . i.e. Eq. (13.71) 
reproduces the even moments = 2, . . . , 14 computed in Refs. Il24ll25ll26l . while Eqs. (13.81) and 
(13.91 ) reproduce the odd moments N = 1, . . . , 13 also obtained in Ref. ll26l . 




5 10 15 5 10 15 

N N 



Figure 1: The perturbative expansion of the anomalous dimension Yns(A^) for four flavours at 
tts = 0.2. In the right part the leading A^^-dependence for large A^^ has been divided out, and the 
corresponding asymptotic limits are indicated as discussed in the text. 



The results (13. 41) . (13.51) and (13.71) for J^siN) are assembled in Fig. 1 for four active flavours 
and a typical value tts = 0.2 for the strong coupling constant (recall that the terms up to order 
a""*"^ are included at N"LO). Numerically, the colour factors take the values Cp = 4/3, = 3 and 
d'^^^dabc/nc = 40/9. Note that the latter normalization is different from that employed in Ref. ll5Sll . 

The NNLO corrections are rather small under these circumstances, amounting to less than 2% 
for N >2. hi large the anomalous dimensions behave as 

y^^^^N) = A„(lnA^ + Y,)-5„-C„i^ + (^i^ (3.10) 

where is the Euler-Mascheroni constant and the coefficients are specified in the next paragraph. 
Thus = y+/lnA^, also shown in Fig. 1, approaches a constant for N ^ oo. The asymptotic 
results are indicated by replacing yns^^(A'^ = 15) by yl^J^{N °°) for the respective highest term 
included in the curves (e.g., for n = 2 at NNLO). Obviously the approach to the asymptotic limit is 



10 



very slow. Yet the results at A'^ — > 0°, which can be derived much easier than the full A/^-dependence 
ll73l . do provide a reasonable first estimate of the corrections. 

The leading large-A'^ coefficients A„, which are also relevant for the soft-gluon (threshold) re- 
summation |l34ll35l|36l|33, are given by 

Ai = ACf 

'67 



A2 = 8C, 



F 



^3 = 16C.Ci(f -^C. + H^3^^^,^)^16C^.^(_| + 2C3) 

+ l6CFCAnf(^-^ + ^^2-\^^+ieCFnj(^-^^ . (3.11) 

The Hy^-independent contribution to the three-loop coefficient A3 is also a new result of the present 
article. Inserting the numerical values of the ^-function and the QCD colour factors it reads 



A2,\nj=Q — 1174.898, in agreement with the previous numerical estimate of Ref. f3T|. The con- 
stants Bn can be read off directly from the terms with 5(1 —x) in Eqs. (14.51) . (14.61) and (14.91) below. 
Surprisingly, the coefficients C„ in Eq. (13.101) . which are also best determined using those x-space 
results, turn out to be related to the A„ by 

Ci = , C2 = 4CfAi , C3 = 8CfA2 . (3.12) 

Especially the relation for C3 is very suggestive and seems to call for a structural explanation. 



4 Results in x-space 

The splitting functions P^s ' (x) are obtained from the N-space results of the previous section by 
an inverse Mellin transformation, which expresses these functions in terms of harmonic polylog- 
arithms Il62ll63ll64l . The inverse Mellin transformation exploits an isomorphism between the set 
of harmonic sums for even or odd and the set of harmonic poly logarithms. Hence it can be per- 
formed by a completely algebraic procedure ll45ll64ll . based on the fact that harmonic sums occur 
as coefficients of the Taylor expansion of harmonic poly logarithms. 

Our notation for the harmonic polylogarithms //mi,...,m„,('^), fnj = 0, ±1 follows Ref. ll64l to 
which the reader is referred for a detailed discussion. The lowest- weight (w = 1) functions H,n{x) 
are given by 

Hq{x) = \nx , H±i{x) = T ln(l Tx) . (4.1) 
The higher- weight (w > 2) functions are recursively defined as 

{■^In'^x, if m^,...,m^ = 0,...,0 

(4-2) 
dzfmi{z)H,n2,...,mM) ' ^IsC 
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with 



1 

X 



/±iW 



1 

1 =Fx 



A useful short-hand notation is 

^o,...,o,±i,o,...,o,±i,...W = ^±(m+i),±(«+i),...W 



(4.3) 



(4.4) 



For w < 3 the harmonic poly logarithms can be expressed in terms of standard poly logarithms; a 
complete list can be found in appendix A of Ref. lH51l . All harmonic polylogarithms of weight 
w = 4 in this article can be expressed in terms of standard polylogarithms, Nielsen functions [74l 
or, by means of the defining relation (14. 2I) . as one-dimensional integrals over these functions. A 
Fortran program for the functions up to weight w = 4 has been provided in Ref. ll75l . 



For completeness we recall the one- and two-loop non-singlet splitting functions 



pS\x) = Cf(27?qq(x)+35(1-X) 



(4.5) 



and 



67 1 1 

P^l^^ix) = AC^Cp (pqq(x) - C2 + yHo + Ho,o 



+ Pqq{-X) 



C2 + 2H_i,o-Ho,o 



+ y(l-x)+5(l-x)[^ + y^2-3C3 )-4Cpnf(pqq{x) 



9 3 " 



+ 3(1-) 



+ 5(l-x) l + ^^2])+4C/(2;.qq(x) 



Hi,o-^Ho-fH2 



-Hon 





r 3 1 




-{1-x) 


_1--H„_ 


-Ho 



2j!7qq(-Jc) ^2 + 2H_l,0 

(4.6) 



pI^^Ux) 



+ l6Cp [Cp-^){ Pqq(-X) ^2 + 2H_i,o - Ho,o 



-2fl-jc) 



-(l+.^)Hoj. (4.7) 
Here and in Eqs. (14.91 1 - (14.1 II) we suppress the argument x of the polylogarithms and use 

Pqqix) = 2{l-x)-'-l-x. (4.8) 

All divergences for x ^ 1 are understood in the sense of +- distributions. 

The three-loop splitting function for the evolution of the 'plus' combinations of quark densities 
in Eq. (12.21 ). corresponding to the anomalous dimension (13.81 ) reads 

Pi?+(x) = 16QC^«/^77qq(x)[^^2-^-9C3-^Ho + 2HoC2-7Ho,o-2Ho,o,o 



1 



g^44V / L 3 

3. 5, 



36 



18 



+ 3Hl,o,0-H3 +-Pqq(-X) -C3-tC2-H_2,0-2H_i^2- 



10, 



-H 



-1,0 



H 



-1,0,0 



1^5 1 , 257 43 1 

+ 2H_i,2 + -HoC2 + ^Ho,o + Ho,o,o-H3 +{l-x) 7^2 - ^ - t^Ho - -Ho,o -Hi 
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1, 1 c.. xfS 167 

• z . 3-C2 + Ho + -Ho,o + 6(l-.)[--- 

+ lec^C/ (pqq W [^^3 - ^^2' - H_3,o - 3H_2C2 - 14H_2,-1,0 + 3H_2,0 + 5H_2,0,0 
151„ 41„ ^ 17„ ^ 13 ^ " 



-5 167^ 1 2 25^ 
.4-^^2 + 20^2 +18^3 



151 41 17 13 23 2 

^^Ho + Y^Ho^2 - ^^0^3 - ^Ho,o - 4Ho,oC2 - y^Ho,o,o + 5Ho,o,o,o + 

+ 8Hij^o,o — 8Hi,3 +H4 



4H_,,-— H„ + - 

67 ^ ^ 31 
+ "^Hi^o — 2Hi^oQ2 + -^"Hlo.o + 1 lHi,o,o,o 

, \ r 1 i- 2 67 31 

+ Pmi-^) [4^2^ - -g^2 + -^^3 + 5H_3,0 



-24HiC3-16Hi,_2,o + y 

67 11 
+ y H2 - 2H2C2 + y H2,o + 5H2,o,o + H3,o 

-32H_2C2-4H_2,-i,o-^H_2,o + 21H_2,o,o + 30H_2,2-^H_i^2-42H_,^, , . 

134 

— 4H_i _2,o + 56H_i _i^2 — 36H_i _i^o,o — 56H_i _i^2 ^H_i^o — 42H_i^oC2 — H3^o 

31 31 13 29 67 

+ 32H_i,3 — — H_i^o,o + 17H_i.o,o,o + ^H_i 2 + 2H_i,2,o + t7:HoC2 + ttHqCs + 77^0,0 

89 31 1 rl33 

+ 13Ho,oC2 + -prHo,o.o-5Ho,o,o,o-7H2C2-^H3-10H4 +(l-x) +4Ho,o.o,o 

167 77 209 
-4—^3 — 2130^3 — 2H_3^o + H-2C2 + 2H_2 -1,0 — 3H_2,o,o + yHo,o,o 5"^^ —7131^2 

31H_i^2-14H_i,_i,o-yH_i,o 



14 

+ 4Hi,o,o + y 



+ VHi,o 



25, 



+ 24H_ 1,2 + 23H_ 1,0,0 + ^ 

1 1 2 

+ 2H2,o,o — 3H4 —5^2—2^2 + 



-43 9 

yC2-3C22 + yH_2,0 

1457 1025 155 " 

tin hlo.O :^n2 + rl2'32 

D 



— HnC? -|- 5Hn nC? H Hn 

2 u,us2T u 26 

-7H_2,o 



-I5H3 



242 



Ho 



-Ho^3 - yHo^2 g- 

151 „ ^ 205^ 
^ + ^2^3 -^^2 



247^2 2ir 15 
--60"^2 +7^C3 + y 



■^5 



3 

+ H-3,o + 4H_2,-i,o — 2^-2,0 — H_2,o,o 

) + 9Hi^3 + 6H1 _2,0 — Hi^oC2 — 



- Ho,o,o,o - 

+ |iH3+H4l +Pqq{-X) ^C2-C2'-^C3-H_3.0 + 8H_2C2 + ^H_2,0-4H_2,0,0 

1 1 67 

- 3H_ 1 ,0,0,0 + y H_ 1 ^2 + 1 2H_ 1 ^3 - 1 6H_ 1 1 ^2 + 8H_ 1 1 ,0,0 + 1 6H_ 1 1 ,2 + y H_ 1 ,0 

11 11 11 67 



2 ■ 50C3-2H_3,o 

, 28„ , ^„ , , rl51 , „ . 

24 

^ 104 

31tt .TT 389 
-— HoC2 + 4HoC3 + ^ 

31 



185 28 

-^Ho,o-22Ho,o,o-4Ho,o,o,o + yH2 + 6H3+5(l -..^ 1^ . 2 

r245 67^ 12^ 2 1. 1043 
^-T8^2 + yC2 +2^3 + ^ 



2Ho,oC2 + —7- 

D 

+ 16 Q^^^^ (pqq W - Yg 

^"2,0,0 



Ho 



-Ho,o - 



■H-2,0 — 4H_2,o,o 



-j-^C2 - ^2^ - y ^3 - H_3,o + 8H_2C2 + y 



-8H_2,2 + 11H 

31 11 
— 3Ho,oC2 — T7:Ho.o,o + Ho,o,o,o + 2H2^2 + -^H3 + 2H4 

1 1^1 523 . 13 5 

2H-3,o - 2^-2^2 + 2^-2,0,0 + ^^"^0 + ^0^3 - y ^0,0 - 2 

— 5H_i,o,o 



y 

11 11 3 1 67 

-1,0^2 + -^H_i,o,o — ^H_i,2 — 8H_i,3 — -Ho — 7Ho^2 — 4Ho^3 — 73-Ho,o 
D i 4 lo 

11 1 r 1 883 1 

Ho.o,o + Ho,o,o,o + 2H2^2 + -^H3 + 2H4 +{l—x) — — — -Hq 0.0,0 + HHi 

J L lUo z 



1 1 

- H_2,-i,o + 2^-3,0 - 2 

8H_iC2 + 4H_i,_i,o + ^H_i,o 



-2Hi,o,o +(1+^: 



' L IU8 2 ' ' 

13 5 

y Ho,o — 2^''^o,o + 2Hi ^2 

-6H_i,2-yC2 + ^C22 



13 



43^ 5 n I 5 ^ 1 ^ 3 1 1 „ 

- - 2^-2fi " ^^0^2 - - ^Ho,os2 + 7H2 - ^H2,o,o + 3H3 + -H4 + -Ho,o^2 

+ ^^2^ - 2 ^2 + + H-2,0 - + - H0C3 + ^^0,0 + 2^0,0,0 + 2^0,0,0,0 

rl657 28U U2 97^ 5^ ]\ 2/ 1 . x Ttt 1 5^ ■ 

I - + 8 C2 + y ^3 - 2 ^5] ) + l6C,n} {-p,,i.) [Ho, - 3 + 3H0 



-^(^-"^im;-^ o .... 

■ - + ^^3] ) + 16Q\ (^Pqq W [5^3 - 4Hl 

1 9 r*; 



+(1- 



rl3 1 



8 

2/1 5.. 1, 



'L54 9 "J ' 'L144 27^" 9^"^J/ ' jyy^^^'y^-^ 

55 5 ^3 10 10 1 2 , ,r5, 

— + -Ho + H0C2 + 2Ho,o - Ho,o,o - y Hi,o - y H2 - 2H2,o - 2H3J + -Pqq(-^) [-^2 

2^3 +H-2,o + 2H_i^2 + yH- 1,0 + H_ 1,0,0 — 2H_i^2 — 2^oC2 — ^Ho.o — Ho,o,o + H3 

rin 1Q 4 9 A 1 vA 95 1 1 



2 

-(1-JC 



+ ^Ho,o + ^H2-5(l-. 



■4 25 1 1 2 

:;-H-i,o - ^Ho + 2^0'O'Oj + 9^0 



[^ + y3Ho,„-lH, + |H,,„ + lH.]+(l+.)[l 

■23 5 , 29, 2 17, n ,/ , ,r 9 >• 



9 2 
lTo^^ -2H_3,o 

+ 6H_2C2 + 12H_2,-i,o - 6H_2,o,o - - ^HoC2 + H0C3 + ^Ho.o - 2Ho,o.o,o + 8H1.3 

Id z 

+ 12Hi^3 + 8Hi,_2,o - 6Hi,o,o - 4Hi,o,o,o + 4Hl2,o - 3H2,o + 2H2,o,o + 4H2,lo + 4H2,2 

' [^C2^ - ^^3 - 6H_3,o + 32H_2C2 + 8H_2,-i,o + 3H_2,o 



J. X , J. 

+ 4H3,0 + 4H3,1+2H4 +77qq( 

— 26H_2,o,o 
+ 48H_i,_i, 



,2 + 40H_i,oC2 + 3H_i,o,o 



2 

1,-1^2 + 40H_i_i,o,o 
-6H_i,2-4H_i,2,o 



- 28H_2,2 + 6H_i^2 + 36H_i^3 + 8H_i,_2,o 

-32H_i,3-^Ho 



— 22H_ 1,0,0,0- 

^HoC2 - I3H0C3 - 14Ho,oC2 - ^Ho,o,o + 6Ho,o,o,o + 6H2C2 + 3H3 +2H3,o + I2H4 

r 31 5 1 

2H_3,o - y + 4H_2,o,o + Ho,oC2 - 3Ho,o,o,o + 35Hi + 6H1 ^2 - Hi,o + -H2,o 

r37 93 81 

-^^2^ - y C2 - y ^3 - 15H_2,o + 30H_iC2 + 12H_i,_i,o - 2H_i,o - 26H_i,o,o 

+ ^H2 - 3H2,o,o - 2H3,o + I3H3 



+ (i-x; 

+ (l+x 



-24H_i,2- — 
Id 

-H4 



539 , 191 8!) 

-— Ho-28HoC2 + ^Ho,o + 20Ho,o,o + -r-. 
Id o 4 

y Ho + 6H0C3 + I9H0C2 - 25Ho,o - 17Ho,o,o 
9.. 18,. o 17,. „ i\ 



+ 4^2 + 33^3 + 4H_3,o 

9. 18..,17^,_„^^ 



D/ 

+ 10H_2,o + yHo 

r29 9 18 17 

- 2H2 - H2,o - 4H3 + S(l -^) - 2C2C3 +3^2 + y C2' + y 



(4.9) 



, ---- 5 4 - "-jy 

The jc-space counterpart of Eq. (13.81) for the evolution of the 'minus' combinations (12.21) is given by 

= /'i?+(x) + 16QC^(c^-^)(;?qq(-x) 

/I A 
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■C2-4C2'-llC3-4H_3,o 



22 44 

+ 32H_2C2 + yH-2,0 - 16H_2,o,o - 32H_2,2 + yH_i^2 + 48H_iC3 - 64H_i,_iC2 

OAQ 22 4^ 

1-1,0^2 + yH-1,0,0 - 12H_i„, - — 

11 



268 22 44 

+ 32H_i,_i,o,o + 64H_i,_i,2 + y-H-i,o + 44H_i,C2 + yH-1,0,0 - 12H_i,o,o,o - yH_i, 

2 134 31 

- 32H_i,3 - 3Ho - -H0C2 - I6H0C3 - y-Ho,o - 12Ho,oC2 - y Ho,o,o + 4Ho,o,o,o + 8H2C2 

14 



22 

+ _H3 + 8H4 



367 1 

+ + 2H_3,o — 2H_2C2 — 4H_2,-i,o — 10H_2,o — 2Ho,o 
140, 



+ 2H_2A0 + 2HoC3 +Ho,oC2 - Ho,o,o,o + 8Hi^2 + —Hi J +{l+x) [32H_i^2 - 18^2 
- 23^3 + y H_i,o - 16H_i,o,o - 32H_i,2 - ^Hq - 29H0C2 + 5Ho,o,o + 24H3 + y H2 



■ 2^2 - 2^3 + 32Ho + I4H0C2 + 2Ho,o,o - I6H3 j + leCpUf [C^ 



Pqqi-X) 2^3 



4 



+ {\-x) 



61 8, 



-H, 



8 41 4 

2Ho,o — 2^-1,0 + y Ho — -H2 



+ 16C^2 (^Q - ^ j [pc,q{-x) [9^3 - 7^2' + 12H_3,o - 64H_2C2 - 16H_2,-i,o - 6H_2,o 

+ 52H_2A0 + 56H_2,2 - 12H_i^2 - 72H_iC3 - 16H_i,_2,o + 96H_i,_i^2 - 80H_i,_ 1,0,0 

- 96H_i,_i,2 - 80H_i,oC2 - 6H_i,o,o + 44H_i,o,o,o + 12H_i,2 + 8H_i,2,o + 64H_i,3 + 3Ho 

+ 3H0C2 + 26H0C3 + 28Ho,oC2 + 9Ho,o,o - 12Ho,o,o,o - I2H2C2 - 6H3 - 4H3,o - 24H4 

- ( 1 - [15 + 8H_3,o + 8H_2,o,o + 6IH0 + 6H0C3 + 2Ho,oC2 - 6Ho,o,o,o + 12Hi ^2 + 6OH1 



+ 8Hi,o 



+ {\+x) 24^2 + 57^3 + 1 OH_2,o - 48H_ 1 ^2 - 4H_ 1 ,0 + 40H_ 1 ,0,0 + 48H_ 1 ,2 



+ 59H0C2 - 22Ho,o - 35Ho,o,o - 22H2 - 4H2,o - 44H3J + 8^2 - 42^3 - 4H_2,o + 42Ho 
-38HoC2 + 14Ho,o-16H2 + 26Ho,o,o + 24H3) . (4.10) 



Finally the Mellin inversion of y^ns^{N) in Eq. (13.91 ) leads to the following result for the leading 
(third-order) difference P^s {x) of the 'valence' and 'minus' splitting functions: 



,(2)s/ 



16n 



! V2 HT 12 ~ 4 ~ ''^ H_24,2 - H_2,o,0 + ^H3 

3 13 91 1 1 r 3 3 

+ 2H_2,-i,o + 2Ho,oC2 - 2H1C2 - ^Hi,o,o + y^HiJ + +^) [H-i,-i,o - 2H-1C2 + ^Ho 

13 1 3 9 29 41 1 

- -g-H_i,o + 2H-i,o,o + 2H_i,2 - -H_2,o + ^HoC2 + ^^^0,0 + -^^2 - H2C2 - 2H2,o,o 



3 

+ 2H4 



1 /I 



3 \x 



+ x' 



3H_ 1 ^2 + 2H_ 1 1 ,0 - 2H_ 1 ,0,0 - 2H_ 1 ,2 + H 1 C2 

2 



1 2 



5C3-2H3 



1 91 9 

+ 2H_2,o + 4H0C2 - 2Ho,o,o + 2Hi ^2] + 24^0 + ^3 - 2 C2 + C2' - H0C3 - H0C2 - 2Ho,oC2 



3 11 

+ g Ho,o - ;^Ho,o,o + 2Ho,o,o,o + H-2,0 - H3 ] . 



(4.11) 



Of particular interest is the end-point behaviour of the harmonic polylogarithms at x — > or 
X ^ I, where logarithmic singularities occur. In the limit x — > 0, the factors \nx are related to 
trailing zeroes in the index field, whereas in the limit x ^ I factors of ln(l — .x:) emerge from 
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leading indices of value 1 . In both limits, the logarithms can be factored out by repeated use of the 
product identity for harmonic polylogarithms, 

E 



Hr {x) 



(4.12) 



Here l+l H,, represents all mergers of fh^ and H,, in which the relative orders of the elements of 
m^v and Hy are preserved. All algorithms for this algebraic procedure have been coded in FORM, 
some explicit examples are given in Refs. lIMlIT^ . 

The large-x behaviour of splitting functions P^}^'^ {£) reflects the large-A'^ behaviour of the 
corresponding anomalous dimensions in Eq. (13.101) . Specifically, the (identical) large-;c behaviour 
of P^i}^'^ [x] is given by 



p(2)±,v 



A3 



+ 538(1 -x) +C3ln(l-x) + 0(1) . 



(4.13) 



The constants A3 and C3 have been specified in Eqs. (13.1 II) and (13.121) . respectively, while the 
coefficients of 5(1 —x) are explicit in Eq. ( 14.91) . At small x the splitting functions can be expanded 
in powers of Inx. For the three-loop non-singlet splitting functions Pn's^^'^l-^) one finds 



P^^^ix) = Dl^Wx + D\\n^x + D^ln^x + D^lnx + 0(1) . 



(4.14) 



Generally, terms up to In^^x occur at order ocf^^ Keeping only the highest n+1 of these, one 
arrives at the N"Lx small-x approximation. Like the large-x coefficients, these contributions can 
be readily extracted from our full results using Eq. (14.121 ). For Pjls^^ we obtain 
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-^Cpflf^-CpHf 
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CpCj^n 



88^2 88^ 
— Cpnj:+—Cp'n 



27 



^pnf 



or, after inserting Ca = 3 and Cp = 4/3 and the numerical values of ^2 and ^3, 
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Dt = 
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295.042 - 32. 1975 rif + 0.592592 n 
1261.11- 152.597/7^ + 4.345679/1 



(4.16) 
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The corresponding coefficients for ' are given by 
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Uq — 

D]^ = 

D2 = 

D — 

^3 — 



1.43210 

35.5556 -3.1 6049 
399.205 - 39.7037 Uf + 0.592592 nf 
1465.93- 172.693 ny + 4.345679 «^ . 



(4.18) 



The coefficients Dq of the leading logarithms in Eqs. ( 14.151) and (14.171 ) agree with the predictions 
in ref. II5T1 based of the resummation of Ref. fj^. Finally the small-jc expansion of Pas^^{x) reads 



^0 



abc 



n 

d^^'d 



-/3 



abc 



n 

d^'^'d 



abc 



n 

d^^'d 



abc 



rir 



^^(18-10^2) 
^^(56 + 2^2-16^3: 



(4.19) 



or, inserting the QCD value of 40/9 for the group factor d"^^'dabc/>^c. 



^ +1.48148 Hy , 
D| ^ +6.89182/1^ , 



Dl = -2.96296 
^ + 178.030 



(4.20) 



The and parts of the functions Pns {^) in Eqs. (14.91) and (14.101) are separately shown 
in Figs. 2-4 together with the approximate expressions derived in Ref. mainly from the 
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integer-A'^ results of Refs. lEHESlEEIl- Also shown for the non-fermionic contributions in Figs. 2 
and 3 are the successive approximations by small-x logarithms as defined in Eq. (14.141) and the text 
below it. As can be seen from Eqs. (14.161) and (14.181) . the coefficients of In^x for Pns^^ increase 
sharply with decreasing power k. Consequently the shapes of the full results in Figs. 2 and 3 are re- 
produced only after all logarithmically enhanced terms have been included. Even then the small-x 
approximations underestimate the complete results by factors as large as 2.7 and 2.0, respectively, 
for Pils^^ and Pns^ atx = 10^^, rendering the small-x expansion (14.141) ineffective for any practi- 
cally relevant value of x. Keeping only the Lx (In'^x) or NLx (In^x and In^x) contributions leads 
to a reasonable description only at extremely small values of x. Therefore, meaningful estimates of 
higher-order effects based on resumming leading (and subleading) logarithms in the small-;c limit 
appear to be difficult. 

The new three-loop nj contribution P^s^^ with the colour structure d'^'^^dabc/nc is graphically 
displayed in Fig. 5 for «/ = 1 . Rather unexpectedly, also this function behaves like In^x for x^Q, 
and here the leading small-x terms do indeed provide a reasonable approximation. In fact, this 
function dominates the small-x behaviour of the non-singlet splitting functions, for rif = A being, 
for example, about 7 times larger than p'^^^{x) aXx= \Q^^ . The presence of a (dominant) leading 
small-x logarithm in a term unpredictable from lower-order structures appears to call into question 
the very concept of the small-x resummation of the double logarithms CX^+^ In^^x. 

In view of the length and complexity of the exact expressions for the functions P^'^\x), it 
is useful to have at ones disposal also compact approximate representations involving, besides 
powers of x, only simple functions like the +-distribution and the end-point logarithms 

©o = l/(l-x)+, Li=ln(l-;c), Lq = Injc . (4.21) 

(2) + 

Inserting the numerical values of the QCD colour factors, P^s in Eq. (14.91) can be represented by 

p'i^^{x) = +1174.898 ©0+1295.3845(1 -x) +714.1 Li + 1641.1 - 3135;c + 243.6x2 

- 522. 1 x^ + LoLi [563.9 + 256.8 Lq] + 1258 Lq + 294.9 + 800/27 + 128/8 1 
+ (^- 183.187 2)0- 173.927 5(1 -x) - 5120/81 Li - 197.0 + 381. Ix + 72.94x2 

+ 44.79 - 1 .497 xL] - 56.66 LqLi - 1 52.6 Lq - 2608/8 1 - 64/27 ) 

+ (-©o-(51/16 + 3C3-5C2)5(1-jc)+jc(1-jc)-1Lo(3/2Lo + 5) + 1 

+ (l-x)(6+ll/2Lo + 3/4L^)) 64/81 . (4.22) 

/j\ ^ 

A corresponding parametrization of P^s in Eq. (14.101) is given by 

Pi^^^{x) ^ + 1 174.898 'Do + 1295.470 5( 1 - x) + 7 14. 1 Li + 1 860.2 - 3505 x + 291. Ox^ 

- 433.2 x^ + LqLi [684 + 25 1 .2 Lq] + 1465 .2 Lq + 399.2L^ + 320/9 + 1 1 6/8 1 

+ f- 183. 187 ©0- 173.933 5(1 -x) -5120/81 Li -216.62 + 406.5 jc + 77. 89 jc^ 
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Figure 2: The -independent three-loop contribution q{x) to the splitting function Pj^{x), mul- 
tiplied by (1 — for display purposes. Also shown in the left part is the uncertainty band derived 
in Ref. ll29ll from the lowest six even-integer moments ll24l l25l l26ll . In the right part our exact 
result is compared to the small-x approximations defined in Eq. (14.141) and the text below it. 




Figure 3: As Fig. 2, but for the splitting function P„g (x). The first seven odd moments underlying 
the previous approximations ll29ll also shown in the left part have been computed in Ref. ll26ll . 
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the uncertainty bands of Ref. ll^ based on the integer moments calculated in Refs. lO^ESl I 
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Figure 5: The first non-vanishing contribution P^^ (x) to the splitting functions P^^ix), compared 
to the approximations of Ref. [291 (where, assuming the completeness of the resummation 0U1I3T1I . 
the possibility of a In^x term was disregarded) and to the small-x expansion in powers of ln.x;. 
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+ 34.76x3 - 1. 136xL^ - 65.43 Lq^i - 172.69 Lq - 3216/81 - 256/81 j 
+ n/(-!Do -(51/16 + 3^3 -5^2) 5(1 -jc)+jc(1-x)-1Lo(3/2Lo + 5) + 1 

+ (l-x)(6+ll/2Lo + 3/4L^)) 64/81 . (4.23) 

Finally the splitting function pji^^ in Eq. (14.1 II) can be approximated by 

Pi^'>\x) ^ Hf (^[Li(-163.9jc"^ -7.208x) + 151.49 + 44.51 jc-43.12x2 + 4.82jc3] [1 -x] 

+ LoLi[-173.1+46.18Lo] + 178.04Lo + 6.892L^ + 40/27[L^-2L^]) . (4.24) 

The identical nj parts of P^^^, the +-distribution contributions (up to a numerical truncation of 
the coefficients involving ^/), and the rational coefficients of the (sub-)leading regular end-point 
terms are exact in Eqs. (14.221) - (14.241) . The remaining coefficients have been determined by fits to 
the exact results, for which we have used the FORTRAN package of Ref. ll75l . Except for x values 
very close to zeros of pji^'{x), the above parametrizations deviate from the exact expressions by 
less than one part in thousand, which should be sufficiently accurate for foreseeable numerical 
applications. For a maximal accuracy for the convolutions with the quark densities, also the co- 
efficients of 5(1 — jc) have been slightly adjusted, by 0.02% or less, using low integer moments. 
Also the complex-A^ moments of the splitting functions can be readily obtained to a perfectly suffi- 
cient accuracy using Eqs. (14.221) - (14.241) . The Mellin transform of these parametrizations involve 
only simple harmonic sums Sm>o{N) (see, e.g, the appendix of Ref. ll60l ) of which the analytic 
continuations in terms of logaritmic derivatives of Euler's F-function are well known. 



5 Numerical implications 

In this section we illustrate the effect of our new three-loop splitting functions /'ns^^'^(jc) on the 
evolution (12.61) of the non-singlet combinations q^s^ix^/jj) of the quark and antiquark distributions. 
For all figures we employ the same schematic, but characteristic model distribution, 

xqt^'{x,4) =x%l-x)'' (5.1) 

with 

a = 0.5 , b = 3 , (5.2) 

facilitating a direct comparison of the various splitting functions contributing to Eq. (12.61) . For the 
same reason the reference scale is specified by an order-independent value for the strong coupling 
constant usually chosen as 

a,i4) = 0.2 . (5.3) 

This value corresponds to /Jq ^ 25. ..50 GeV^ for as(M|) = 0.114. ..0.120 beyond the leading 
order, a scale region relevant for deep-inelastic scattering both at fixed-target experiments and. 
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for much smaller x, at the ep collider HERA. Our default for the number of effectively massless 
flavours is n^: = 4. The normalization of is irrelevant for our purposes, as we consider only the 
logaritmic derivatives 4ns = dlnq^^/dln/jj. 




Figure 6: The perturbative expansion of the logarithmic scale derivative dlnq:^^/dln^j for a char- 
acteristic non-singlet quark distribution ;c^+ = x^'^{i —x)^ at the standard scale /Ur = l^f- 

The scale derivatives of the three non-singlet distributions are graphically displayed in Figs. 6 
and 7 over a wide region of x. At large x the NNLO corrections are very similar in all cases, 
amounting to 2% or less for x > 0.2, thus being smaller than the NLO corrections by a factor of 
about eight. The same suppression factor is also found for ^ns(-^) the region 10^^ ^ x < 10^^. 
The NNLO effects are even smaller for ^+ at small x, but considerably larger for at x < 10^-^. 
For example, at x ~ 10^"*, where P^^^ {x) exceeds p'^ {x) by a factor of about 8 as discussed in 
the paragraph above Eq. (14.211 ). the ratio of the corresponding corrections in Fig. 7 amounts to 2.5. 
Recall that the scale derivatives ( 12.61 ) do not probe the splitting functions locally in x due to the 
presence of the Mellin convolution. 

The numerical values for q^^{x,^Q) are presented in Tab. 1 for four characteristic values of x. 
Also illustrated in this table is the dependence of the results on the shape of the initial distribution, 
the number of flavours and the value of the strong coupling constant. The relative corrections 
are rather weakly dependent of the large-jc power b in Eq. (15.11) . They increase at small x with 
increasing small-x power a, i.e., with decreasing size of q^^. At large x, where the n^d"^'^ dabc/nc 
contribution P^^ is negligible, the NNLO corrections decrease with increasing rij:. At small-x 
this decrease is overcompensated in ql^ by the effect of P^s- Except for very small momentum 
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Figure 7: As Fig. 6, but for the scale derivatives of the two other non-singlet combinations ^ns'^- 

fractions x < 10^^ (where the non-singlet quark densities play a minor role for most important 
observables) the NNLO corrections amount to 15% or less even for a strong coupling constant as 
large as tts = 0.5. Hence the non-singlet evolution at intermediate and large x appears to remain 
perturbative down to very low scales as used in the phenomenological analyses of Refs. iTTTllTSl 
and in non-perturbative studies of the initial distributions like those of Refs. Il38ll39ll40ll4ni42ll43ll . 

Another conventional way to assess the reliability of perturbative calculations is to investi- 
gate the stability of the results under variations of the renormalization scale /Jr- For /Jr ^ /Jf the 
expansion in Eq. (12.61) has to be replaced by 



(5.4) 



,(2),/ 



,(0) 



(1),/ 



In4 + |3^/'i^n4 



where p^t represent the MS expansion coefficients of the p-function of QCD ir79ll80ll8ni82ll . 

In Fig. 8 the consequences of varying over the rather wide range < ^u^ < are 
displayed for at six representative values of x. The scale dependence is considerably reduced 
by including the third-order corrections over the full x-range. At NNLO both the points of fastest 
apparent convergence and the points of minimal ^u^-sensitivity, dq^^/dfir = 0, are rather close to the 
'natural' choice /Jr = /Jf for the renormalization scale. 
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X 


LO 


NLO 


NNLO 


n 




ri/n 


default (Fig. 7) 


10-4 
0.002 
0.25 
0.75 


6.546-10-2 
5.632-10-2 
-5.402-10-2 
-1.949- 10 1 


8.424-10-2 
6.875-10-2 
-6.331-10-2 
-2.189-10-1 


9.163- 10 2 
7.041 - 10-2 
-6.457-10-2 
-2.222-10-1 


0.287 
0.221 
0.172 
0.123 


0.088 
0.024 
0.020 
0.015 


0.31 
0.11 
0.12 
0.12 


a = 0.8 


10-4 
0.002 
0.25 
0.75 


1.660-10 1 
1.249-10-1 
-4.352-10-2 
-1.930-10 1 


2.351-10-1 
1.583-10-1 
-5.171-10-2 
-2.168-10-1 


2.818-10-1 
1.650-10-1 
-5.283-10-2 
-2.200-10-1 


0.417 
0.268 
0.188 
0.123 


0.198 
0.042 
0.022 
0.015 


0.48 
0.16 
0.12 
0.12 


b = 5 


10-4 
0.002 
0.25 
0.75 


6.474-10-2 
5.324-10-2 
-7.835-10-2 
-2.300-10-1 


8.278-10-2 
6.432-10-2 
-9.022-10-2 
-2.580-10-1 


8.917-10-2 
6.546-10-2 
-9.180-10 2 

-2.619-10 1 


0.279 
0.208 
0.151 
0.122 


0.077 
0.018 
0.018 
0.015 


0.28 
0.09 
0.12 
0.12 




10-4 

0.002 

0.25 

0.75 


6.546-10-2 
5.632-10-2 
-5.402-10-2 
-1.949-10-1 


8.480-10-2 
6.942-10-2 
-6.406-10-2 
-2.219-10-1 


9.187-10-2 
7.174-10-2 
-6.588-10-2 
-2.269-10-1 


0.295 
0.233 
0.186 
0.139 


0.083 
0.033 
0.028 
0.023 


0.28 
0.14 
0.15 
0.16 


fif = ^ and tts = 0.5 


10-4 

0.002 

0.25 

0.75 


1.636- 10 1 
1.408-10-1 
-1.350-10-1 
-4.871-10-1 


2.845-10-1 
2.227-10-1 
-1.978-10-1 
-6.563-10-1 


3.949-10-1 
2.589-10-1 
-2.262-10-1 
-7.346-10-1 


0.739 
0.581 
0.465 
0.347 


0.388 
0.163 
0.144 
0.119 


0.53 
0.28 
0.31 
0.34 



Table 1: The LO, NLO and NNLO logarithmic derivatives q'^^ = d\nq^^/d\n/uj at four representa- 
tive values of x, together with the ratios r„ = N"L0/N"-1L0 — 1 for the default input parameters 
specified in the first paragraph of this section and some variations thereof. 
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Figure 8: The dependence of the NLO and NNLO predictions for q^^ = dlnq^^/dlufjy on the 
renormalization scale /Jr for six typical values of x. The initial conditions are as in Fig.|5l 



The relative scale uncertainties of the average results, conventionally estimated by 

max[qi^{x,i4 = lljj..Aijj)]-mm[qi^{x,i4 = l^^^ 

2 1 average [q}^,{x,^if = ^tuj... 4nj)\ \ 

is shown in Fig. 9 for all three cases i = ±,v. These uncertainty estimates amount to 2% or less ex- 
cept for jc < 10 an improvement by more than a factor of three with respect to the corresponding 
NLO results. Taking into account also the apparent convergence of the series in Figs. 6 and 7, it is 
not unreasonable to expect that the effect of the four- loop non- singlet splitting functions — which 
most likely will remain uncalculated for quite some time — will be less than 1% for x > 10^. 
This expectation is consistent with the Pade estimates of pjs^' employed in Ref. ll83ll for the N^LO 
large-;c evolution of the deep-inelastic structure functions F2 and F3. At very small values of jc the 
higher-order corrections will presumably be considerably larger. 
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Figure 9: The renormalization scale uncertainty of the NLO and NNLO predictions for the scale 
derivative of ^^g, z = ±,y, as obtained from the quantity Aq^^ defined in Eq. (I5.5I) . Here and 
in Figs. 6 and 7 the spikes close to x = 0. 1 reflect the sign-change of q^^ and do not constitute 
appreciable absolute corrections and uncertainties. 



6 Summary 



We have calculated the complete third-order contributions to the splitting functions governing the 
evolution of unpolarized non-singlet parton distribution in perturbative QCD. Our calculation is 
performed in Mellin-A'^ space and follows the previous fixed-A'^ computations Il24ll25ll26ll inasmuch 
as we compute the partonic structure functions in deep-inelastic scattering at even or odd A'^ using 
the optical theorem and a dispersion relation as discussed in ll25l . Our calculation, however, is not 
restricted to low fixed values of but provides the complete A'^-dependence from which the x-space 
splitting functions can be obtained by a (by now) standard Mellin inversion. This progress has been 
made possible by an improved understanding of the mathematics of harmonic sums, difference 
equations and harmonic polylogarithms [591 I^H |35|, and the implementation of corresponding 
tools, together with other new features [i53J . in the symbolic manipulation program Form [»52l 
which we have employed to handle the almost prohibitively large intermediate expressions. 

Our results have been presented in both Mellin-A^ and Bjorken-x space, in the latter case we 
have also provided easy-to-use accurate parametrizations. Our results agree with all partial results 
available in the literature, in particular we reproduce the lowest seven even- or odd-integer mo- 
ments computed before ll24l l25l l26l . We also agree with the resummation predictions ll30l |3T1 
for the leading small-x logarithms In^x of the splitting functions P^{x) and P^{x) governing the 
evolution of flavour differences of quark- antiquark sums and differences. However, an unpre- 
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dieted term of the same size is found also for the new d dahc/^c contributions P^^ to the splitting 
function for the total valence distribution. At large x we find that the coefficients of the leading 
integrable term ln(l —x) at order n is proportional to the coefficient of the (only) +-distribution 
1/(1 — x)+ at order n — \ , a result that seems to point to a yet unexplored structure. 

We have investigated the numerical impact of the three-loop (NNLO) contributions on the 
evolution of the various non-singlet densities. The effect of the new contribution P^{x) is very 
small at large x but rises sharply towards x ^ 0, reaching 10% for a standard Regge-inspired ^/x 
initial distributions at jc ~ 10^. At jc > 10^, on the other hand, the perturbative expansion for 
the scale dependences dlnqji^ix, fjj) / dln^j appear to be very well convergent. For = 0.2, for 
example, the NNLO corrections amount to 2% or less for four flavours, a factor of about 8 less 
than the NLO contributions. Also the variation of the renormalization scale leads to effects of 
about ±2% at NNLO in this region of x. Corrections of this size are comparable to the dependence 
of the predictions on the number of quark flavours, rendering a proper treatment of charm effects 
rather important even for large-jc non-singlet quantities, see Refs. ll84l[85l and references therein. 

Form files of our results, and Fortran subroutines of our exact and approximate ;t:-space 
splitting functions can be obtained from the preprint server http : / / arXiv .org by downloading 
the source. Furthermore they are available from the authors upon request. 
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